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We show that the confining property of the one-gluon propagator, in Coulomb gauge, is linked to the unbroken 
realization of a remnant gauge symmetry which exists in this gauge. An order parameter for the remnant gauge 
symmetry is introduced, and its behavior is investigated in a variety of models via numerical simulations. We 
find that the color-Coulomb potential, associated with the gluon propagator, grows linearly with distance both in 
the confined and — surprisingly — in the high-temperature deconfined phase of pure Yang-Mills theory. We also 
find a remnant symmetry-breaking transition in SU(2) gauge-Higgs theory which completely isolates the Higgs 
from the (pseudo)confinement region of the phase diagram. This transition exists despite the absence, pointed 
out long ago by Fradkin and Shenker, of a genuine thermodynamic phase transition separating the two regions. 



I. INTRODUCTION 

It is well known that the phase transition from a confined 
to an unconfined phase in a non-abelian gauge theory is asso- 
ciated with the breaking of global center symmetry. If a non- 
trivial center symmetry exists and is unbroken, then Polyakov 
line expectation values are zero, and in consequence the free 
energy of a static color charge is infinite in an infinite vol- 
ume. But the confined phase may be associated with other 
global symmetries as well. In this article we will focus on 
the remnant gauge symmetry which is found after imposing 
Coulomb gauge. What is interesting about this symmetry is 
that its unbroken realization implies the existence of a confin- 
ing color-Coulomb potential, and this in turn is a necessary 
(but not sufficient) condition for confinement. 

The color-Coulomb potential arises from the energy of the 
longitudinal color-electric field in Coulomb gauge, and corre- 
sponds diagrammatically to instantaneous dressed one-gluon 
exchange between static sources. We have previously studied 
this potential numerically in pure lattice Yang-Mills theory at 
zero temperature, and found that the potential rises linearly 
with color charge separation [1,2], albeit with a string tension 
'^coui which is significantly higher (by about a factor of three 
[2]) than the string tension o of the static quark potential. In 
this article we introduce an order parameter for remnant gauge 
symmetry breaking and study its behavior, as well as the be- 
havior of the color-Coulomb potential, in 

• SU(2) gauge-Higgs theory, with the Higgs field in the 
adjoint representation; 

• pure SU(2) gauge theory, in the confined and the high- 
temperature deconfined phase; 
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• SU(2) gauge-Higgs theory, with the Higgs field in the 
fundamental representation; 

• compact QED4. 

In the first, second, and fourth theories the action has a non- 
trivial global center symmetry, which may be broken spon- 
taneously in some range of couplings or temperature. In 
the third case center symmetry is explicitly broken, and the 
asymptotic string tension vanishes at all couplings. By study- 
ing these different cases, we can explore to what extent rem- 
nant symmetry breaking is correlated with center symmetry 
breaking, and whether the confining Coulomb potential is al- 
ways associated with a confining static quark potential. 

This article is organized as follows: In section II, below, we 
introduce the remnant symmetry order parameter, and relate 
it to the confining properties of the color-Coulomb potential. 
We also examine scaling of Ocoui with p. Section III concerns 
the gauge-adjoint Higgs model, where we find perfect cor- 
respondence between remnant symmetry and center symme- 
try breaking. But this correspondence is lost already in pure 
Yang-Mills theory at high temperature, studied in section IV, 
where we find that Coulomb confinement and unbroken rem- 
nant symmetry persist in the deconfined phase. A possible 
explanation of this phenomenon is discussed. In section V we 
review the Gribov confinement scenario in Coulomb gauge 
and show that it accords with the vortex dominance scenario, 
by gauge transforming from the maximal center gauge to 
the minimal Coulomb gauge. In section VI we present our 
results for the gauge-fundamental Higgs model, where the 
gauge-Higgs interaction breaks global center symmetry ex- 
plicitly. In this case we find very clear numerical evidence of 
a remnant symmetry-breaking transition which is unaccom- 
panied by a true thermodynamic phase transition, and also 
argue for the existence of such a transition from a lattice 
strong-coupling analysis. This result is in complete accord 
with the earlier work of Langfeld [3, 4], which found rem- 
nant symmetry breaking in Landau gauge, and it implies that 
a sharp distinction can be made between the Higgs and the 
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(pseuclo)confimng regions of the gauge-Higgs coupling plane. 
This distinction exists despite the fact, pointed out by Fradkin 
and Shenker [5], that these regions are continuously connected 
in the usual sense of thermodynamics. The remnant symme- 
try breaking in the adjoint Higgs theory is re-examined in sec- 
tion VII, where we draw some conclusions about the measure 
of abelian configurations in the fundamental modular region. 
Our methods are applied to compact QED4 in section VIII. 

As in the previous work of ref. [2], we also investigate 
numerically the relevance of center vortices to the existence 
of a confining Coulomb potential, particularly in the high- 
temperature deconfined phase of pure Yang-Mills theory, and 
in the (pseudo)confined phase of gauge-fundamental Higgs 
theory. In section V it is shown that thin center vortices lie 
on the Gribov Horizon, which may be relevant to their dra- 
matic effect on the Coulomb potential. Vortex removal, by 
the de Forcrand-D'Elia procedure |6|, is found in every case 
to convert a confining Coulomb potential to a non-confining, 
asymptotically flat potential. 

In section IX, we provide a 2-way translation between the 
temporal gauge (Aq = 0) and the minimal Coulomb gauge. 
This allows our measurements, which are made by gauge- 
fixing to the minimal Coulomb gauge, to be equivalently de- 
scribed in temporal gauge. We show that the state in which our 
measurements are made is, in temporal gauge, a quark-pair 
state of the type introduced by Lavelle and McMullan [7]. It 
has correct gauge-transformation properties, although it does 
not make use of Wilson lines running between the sources. 
Section X contains some concluding remarks. 

II. COULOMB ENERGY AND REMNANT SYMMETRY 

On the lattice, minimal Coulomb gauge consists of fixing to 
the configuration on the gauge orbit maximizing the quantity 

3 

/? = ££Tr[f/;fe(x,0]. (2.1) 

x,« k=l 

Maximizing R does not fix the gauge completely, since there 
is still the freedom to perform time-dependent gauge transfor- 
mations 

Ukix,t) ^ g{t)Uk{x,t)g\t) (A; =1,2,3), 
Uo{x,t) ^ g{t)Uo{x,t)g^t + l). (2.2) 

To understand the role of this remnant gauge symmetry, con- 
sider a state obtained by operating on the lattice Yang- 
Mills vacuum state 'I'o (i.e. the ground state of the Coulomb 
gauge transfer matrix) with a heavy quark operator q" 

^';[x■,A]=q"{x)^'o[A]. (2.3) 

Evolve this state for Euclidean time T <nt, where n, is the lat- 
tice extension in the time direction. Dividing out the vacuum 
factor exp[— £0^] and taking the iimer-product with 1'^[x;A], 
we have 



n 

= re-'"^^{L'"'{x,l,T)) (2.4) 

where index n refers to summation over charged energy eigen- 
states, ^ 00 is the heavy quark mass, F is a numerical fac- 
tor, and L(x.fi .f2) is a Wilson line 

L{x,ti,t2)=Uo{x,ti)Uo{x,h + l)...Uo{x,t2). (2.5) 

In a confining theory, the excitation energy AEn = E„— Eq 
of any state ^„ containing a single static quark must be infi- 
nite. It follows tiiat G'"'{T) = 0, and therefore that flie expec- 
tation value of any timelike Wilson line L(x, 1,7") must vanish. 
Now, under the remnant gauge symmetry, L transforms as 

L{x,ti,t2)-^g{ti )L{x,h ,t2)g\t2 + \). (2.6) 

If r = f2 — f] + 1 is less than the lattice extension in the 
time direction, then not only L but also Tr[L] is non-invariant 
under the remnant gauge symmetry. ' If the remnant symme- 
try is unbroken, then (L) must vanish. The confining phase 
is therefore a phase of unbroken remnant gauge symmetry; 
i.e. unbroken remnant symmetry is a necessary condition for 
confinement. 

There are several issues which require some further com- 
ment. First of all, why is unbroken remnant symmetry not 
also a sufficient condition for confinement? The answer is that 
in an unconfined phase, where there exist finite energy states 
containing a single static charge in an infinite volume, there 
is still the possibility that these finite energy states have van- 
ishing overlap with as defined in eq. (2.3). Thus (L) = 
and unbroken remnant symmetry could be found, in principle, 
also in the absence of confinement. Secondly, it is obviously 
impossible to insert a single charge in a finite volume with, 
e.g., periodic boundary conditions; electric field Unes starting 
from the static charge must end on some other charge, regard- 
less of whether or not the theory is in a confining phase. How, 
then, is this fact reflected in the remnant symmetry-breaking 
criterion? The same question can be raised in connection with 
Polyakov lines, and the answer is the same: Strictly speak- 
ing, spontaneous symmetry breaking cannot occur in a finite 
volume, so (L) = always, consistent with the absence of 
an isolated static charge. Nevertheless (again like Polyakov 
lines), it is possible to construct an order parameter which de- 
tects the infinite-volume transition via finite-volume calcula- 
tions which are subsequently extrapolated to infinite volume.^ 
We will construct such an operator below. Finally, there is the 
question of Elitzur's theorem. Although the remnant symme- 
try is global on a time shce, it is local in the time direction, 
and according to the theorem local symmetries cannot break 
spontaneously. So how could we ever have (L) 7^ 0, even in an 



If r = n,, then L is a Polyakov line, whose trace is invariant under gauge 
transformations, but non-invariant under global center transformations. 
^ In the case of spontaneous center symmetry breaking, the accepted order 
parameter on a finite lattice is the absolute value of the spatial average of 
Polyakov lines. 
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infinite volume? The answer is that in fact the average value 
of L{x,t,T +t) does indeed vanish on an infinite lattice, in ac- 
cordance with the Elitzur theorem, providing the averaging is 
done over all spatial x and all times On a time sUce, how- 
ever, the symmetry is only global, and it is possible in any 
given configuration that the average value of L{x,t,T +t) is 
finite on an infinite lattice, when averaged over all x at fixed 
time t. This is what we will mean by the phrase "spontaneous 
breaking of the remnant symmetry," and it involves no actual 
violation of the Ehtzur theorem. 

With these points in mind, we propose to construct an order 
parameter for remnant symmetry breaking from the timelike 
link variable averaged over spatial volume at fixed time. Let 
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(2.7) 



where = HxHyn^ is the is the 3-volume of a lattice time-shce. 

~ 1 /2 

If remnant symmetry is unbroken, then U (f ) = + (9(1 /V^ ) 
in any thermalized lattice configuration. The order parameter 
Q is defined to be, for SU(2), 
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Q = -L{\/-2nu{t)uf{t)]). 



(2.8) 



Then Q is positive definite on a finite lattice, and on general 
grounds 



Q = c + 



b ^jjgj-g r c = in the symmetric phase, 
JVi \c>Q in the broken phase. 

(2.9) 

If Q extrapolates to a non-zero value as V3 °o, then the rem- 
nant symmetry in Coulomb gauge is spontaneously broken. 

Next, we make the connection between unbroken remnant 
symmetry and the existence of a confining Coulomb poten- 
tial. We first recall that the Hamiltonian operator in Coulomb 
gauge has the form H — Hgi^e + Hcoui where, in the continuum 



glue = \jd^x (i-5£te.«^ .£fr.«j-5 +Bf 



Hcoul = \ j d^xd^yr'^p\x)yK''\x,y-A)9\y)rT- 



K"''{x,y;A) = 



1 



V-Z)(A) 

?abc 



(-V2) 



I 



V-Z)(A) 



ab 



3 = det[-V •Z)(A)] 



(2.10) 



and the factors of J arise from operator ordering consider- 
ations [8]. It is understood that A is identically transverse in 
Coulomb gauge, A=Al^. Note that J commutes with all quan- 
tities except E"" , and in particular with p^(x) and K^^{x.y;A^. 
The expectation-value of K{x^y;A) is the instantaneous piece 
of the (AqAo) gluon propagator, i.e. 



kA%{x)A\{y)) 



D(x-y)5' 



ab 



= D(x-y)5«*5(xo-M) 
-|- non-instantaneous 

1 . 1 



S/-D[A] 



V-£)[A] 



(2.11) 



as shown in [9]. We see from these expressions that 
the Coulomb interaction energy between two charged static 
sources is given by instantaneous (dressed) one-gluon ex- 
change. 

Now consider a physical state in Coulomb gauge containing 
massive quark-antiquark sources 



(2.12) 



which is invariant under the remnant symmetry. The excita- 
tion energy is 



= ^coul{R)+Ese, 



(2.13) 



where Ese is an /?-independent constant, on the order of the in- 
verse lattice spacing, to be specified below. The /^-dependence 
of E can only come from the expectation value of the non- 
local quark-quark part of the Hamiltonian 

Hqq = hj d^xd^y p^(x) K">'{x,y;A) p^(j). (2.14) 

Thus the /^-dependent piece Vcoui{R) can be identified as the 
Coulomb potential due to these static sources. Moreover, the 
same kernel K{x,y;A) appears in Hqq and in the instantaneous 
part D{x) of the (dressed) one-gluon propagator (AqAo). This 
yields the formula, 

VcouM+Ese = C.(z)(0)-D(x)), (2.15) 

where Cr (= 3/4 for the SU(2) gauge group) is the Casimir 
factor in the fundamental representation. 

The correlator of two Wilson lines can be expressed in 
terms of the Hamiltonian operator and the state as fol- 
lows: 



G{R,T) = (iTr[Z,t(x,0,T)L(y,0,T)]) 



(2.16) 



where R = |x — y|, and L is now the timelike Wilson Une in 
the continuum theory. We have 



G{R,T) = '£\{^n\^qq)\^e 



and we define the logarithmic derivative 



V{R,T) = - — log[G{R,T) 



-AE„T 



(2.17) 



(2.18) 



It is easy to see that the Coulomb energy is obtained at T — > 0, 
i.e. 



E = Vcoul{R)+Es, 
= ViR,0). 



(2.19) 



The minimum energy of a state containing two static quark 
antiquark charges, which in a confining theory would be the 



4 



energy of the flux tube ground state, is obtained in the opposite 
T oo hmit 



lim V{R,T) 



(2.20) 



and in this limit V{R) is usual static quark potential. 

The idea of using the correlator of timelike Wilson lines 
to compute the static quark potential, in Coulomb gauge on 
the lattice, was put forward some years ago by Marinari et al. 
[10]. These authors also noted that the remnant symmetry in 
Coulomb gauge is unbroken in the confining phase. 

We now recall an inequality first pointed out by one of us 
(D.Z.) in ref. [11]. With a lattice regularization, Ese and E'^^, 
are finite constants. In a confining theory, both of these con- 
stants are negligible compared to V ) , at sufficiently large R. 
But since E^m < E . it follows that 



V{R)<Vcoui{R) 



(2.21) 



asymptotically. The intriguing implication is that if confine- 
ment exists at all, then it exists already at the level of dressed 
one-gluon exchange in Coulomb gauge. But we also see that 
because the Coulomb potential is only an upper bound on the 
static potential, a confining Coulomb potential is a necessary 
but not a sufficient condition for the existence of a confining 
static quark potential. 

On the lattice, the continuum logarithmic derivative in eq. 
(2.18) is replaced by 



V{R,T) = -log 
a 



G{R,T) 



G{R,T + a) 



(2.22) 



where a is the lattice spacing. In particular, in lattice units 
a = l, 



V{R,0) 



-log[G(/?,l)] 



= -log 



iTr[[/o(x,l)f/J(y,l)] 



(2.23) 



and at large p, where the lattice logarithmic derivative ap- 
proximates the continuum, V{R,0) provides an estimate of the 
Coulomb potential Vcoui{R) (up to an additive constant E^e). 

In eq. (2.23) the relation between the confining property of 
the Coulomb potential, and the unbroken realization of rem- 
nant symmetry, is manifest. For if 2 ^ at infinite volume, 
then also 



lim(Tr[f/o(x,Of/o^y,0]) = 1™ (Tr[[/(Of/^(0]) 



= 



(2.24) 



in which case the potential V{R,0) rises to infinity as /? °o. 
Conversely, if 2 > 0, then the limit in eq. (2.24) is finite, and 
V{R,0) is asymptotically flat. Since Vcoui{R) « V{R,0) is an 
upper bound on the static quark potential, we see again that 
unbroken remnant symmetry is a necessary but not sufficient 
condition for confinement. 

V{R, T) has been computed numerically for pure SU(2) lat- 
tice gauge theory at a range of lattice couplings P in ref. [2]. 
We recall the essential results: 



1. V{R,T) increases linearly with R at large R and all T, 
at any coupling. 

2. The associated string tension o{T) converges (from 
above) to the usual asymptotic string tension a, at any 
given p, as T increases. 

3. At weaker couplings, the Coulomb string tension 
Gcoui = cr(0) is substantially greater (by about a factor 
of three) than the asymptotic string tension. 

4. Removing center vortices from lattice configurations 

(by the de Forcrand-D'Elia procedure [6]) sends 
o{T) ^ at all T, including the Coulomb string ten- 
sion (Seoul ^ at T = 0. 

In the following sections we extend the investigation to 
models including scalar matter fields, and to SU(2) gauge 
theory across the high temperature deconfinement transition.-' 
First, however, we would like to remark on the scaling prop- 
erties of Ocoul ~ (J(0) in pure SU(2) gauge theory at zero tem- 
perature. The ratio of the Coulomb (T = 0) to the asymptotic 
(T 00) string tensions, reported in ref. [2], varies somewhat 
with P in the range (P G [2.2,2.5]) of couplings investigated. 
The ratio a(0)/a tends to rise in this interval, as shown in 
Fig. 1. However, it is known that a does not quite conform 
to the two-loop scaling formula associated with asymptotic 
freedom, in the range of P we have studied, and it is always 
possible that scaling sets in at different values of P for differ- 
ent physical quantities. What we find is that when our values 
for a(0) are divided by the asymptotic freedom expression 



102/121 



exp 



(2.25) 



relevant to the SU(2) string tension, the ratio a(0)/F(P) is 
virtually constant as seen in Fig. 2. This fact suggests that 
scaling according to asymptotic freedom may set in earlier for 
the Coulomb string tension Ocoul ~ <j(0) than for the asymp- 
totic string tension a. 



A. Divergent constant in D{R) 

We would also like to remark, at this stage, on a subtlety 
in identifying the color Coulomb potential with dressed one- 
gluon exchange.^ The energy expectation value of the static 
qq state, E (R) = Vcoui{R) +Ese, is finite for finite quark sepa- 
ration R and finite lattice spacing a; in fact we have calculated 
this quantity (= V{R,Q)) numerically. So it might seem nat- 
ural, from eq. (2.15), to identify the /?-dependent Coulomb 



^ A note on gauge fixing to Coulomb gauge: In this investigation we gener- 
ate eight random gauge copies of each lattice configuration, and carry out 
gauge-fixing on each configuration by over-relaxation for 250 iterations. 
The best copy of eight copies is then chosen, and the over-relaxation pro- 
cedure is continued until the average value of the gauge-fixed links has 
changed, in the last 10 iterations, by less than 2 x 10^^. 
We thank M. Polikarpov for a helpful discussion on this point. 
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string tensions extracted from V(R,0) 



E 




FIG. 1: The ratio o(0) /(J at various P, from fits whicli either include 

(fc = 7r/12) or do not include (b = 0) the Liischer term. This ratio 
should equal the ratio of o^-gui /a in the continuum limit. (Data points 
come from numerical simulations on lattice sizes 16*, 20'*, 20*, and 
32* at p = 2.2,2.3,2.4, and 2.5 respectively.) 
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FIG. 2: The ratio a(0)/F(p) vs. p, again from fits with (b = Jl/12) 
or without (b = 0) the Liischer term. F(P) is the asymptotic freedom 
expression given in the text; a constant ratio implies scaling accord- 
ing to the two-loop beta function. 



subtraction at /? = /?o which removes the infrared divergence, 
i.e. 

VcouliR) = -Cr(D{R)-D{Ro)) 
= ■E{R)-'E{Ro) 

Ese = Cr(D{0)-D{Ro)) 

= £(/?o) (2.26) 

Defined in this way, Vcoui{R) crosses zero at the subtraction 
point, and Ese contains the ultraviolet, but not the infrared, 
contributions to the quark-antiquark self-energies. 

As a check of the cancellation (or non-cancellation) of in- 
frared divergences, consider a colored state consisting, e.g., of 
two static quarks, rather than a quark and antiquark. The en- 
ergy of such a state could be extracted from an LL correlator, 
which is zero if remnant symmetry is unbroken. The energy 
is therefore infinite, and according to our previous analysis 
would be proportional to D{Q) +D{R). In this case, the di- 
vergent constant in D{R) adds to, rather than subtracts from, 
the divergent constant in D{0), and the resulting sum is diver- 
gent, as it should be. The argument can be readily general- 
ized to baryonic states in SU(A^) gauge theories composed of 
static charges. The energy of a color singlet state, with charges 
at points Xi,X2, . . . ,Xn is obtained from the logarithmic time 
derivative of the correlator 

G({x,},r) = e,i..,^e,j...^j,(L"^i(xi,0,r) . . .V^J^xnA t)). 

(2.27) 

The order T contributions to G({x,},r) are terms propor- 
tional to D{0) and D{xm — x„),m^n, with differing signs. On 
the other hand, for a color singlet state, the operator eeL...L 
is a T-independent constant in the xi = X2 = ■ ■ ■ = xn coin- 
cidence limit. From this it is clear that the propagators com- 
pletely cancel in the coincidence limit, and any constant terms 
in the propagators cancel in general. This means that the en- 
ergy of a color singlet baryonic state is finite. Conversely, the 
divergent constants do not cancel in color non-singlet states, 
so their energies are infinite. 



interaction energy as Vcoui{R) = —CrD{R). That cannot be 
quite right, however. The reason is that jCrD{0) is the en- 
ergy of an isolated quark state in an infinite volume; it is the 
energy we would extract from the logarithmic time derivative 
of G(r) = (iTrL(0,0,r)). In the case of unbroken remnant 
symmetry we have G{T) — 0, and therefore ^CrD{0) = oo in 
an infinite volume, even though the lattice spacing a is non- 
zero. This infinity, which has a non-perturbative, infrared ori- 
gin, should not be confused with the usual ultraviolet contri- 
bution to the quark self-energy, which is only infinite in tiie 
continuum limit. Since E (/?) is finite, the infrared divergence 
in CrD{0) must be cancelled, in eq. (2.15), by a correspond- 
ing divergent constant contained in CrD{R). In other words, 
only the difference D{0) - D{R) is finite, and Vcoui{R), if fi- 
nite, differs from —CrD{R) by an infinite constant. In order 
that Va,ui{R) and E^e be separately finite and well-defined, we 
may relate them to the gluon propagator with an (arbitrary) 



ni. SU(2) GAUGE-ADJOINT fflGGS THEORY 

The lattice action for SU(2) gauge theory with a Higgs field 
in the adjoint representation of the gauge group is 

5 = P £ iTr[t/J/f/+f/^] (3.1) 

plaq 

+ lY.r{x)<^''{x+n)Tr[a"U,{x)o'U^^{x)] 

with the radially "frozen", three-component Higgs field (j) sub- 
ject to the restriction 

£rWr(^) = l- (3.2) 

a=l 
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FIG. 3: Phase diagram of the SU(2) adjoint Higgs model. The pla- 
quette energy Ep and the remnant symmetry order parameter Q locate 
the same transition line between the confined and Higgs phases. 
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FIG. 4: Plaquette energy Ep vs. P at three values of y in the gauge- 
adjoint Higgs model. 



This theory was first studied numerically by Brower et al. [12]. 
In addition to SU(2) gauge symmetry, the action is also invari- 
ant under global Z2 center symmetry 



(3.3) 



for some choice oft= to, with all other fields unchanged. The 
existence of this apparently innocent global symmetry in the 
action has profound consequences; in its absence the static 
quark potential is asymptotically flat, and there can be no truly 
confined phase. The significance of center symmetry to the 
confinement property in general is reviewed in ref. [13]. 

Since the action (3.2) is symmetric under global center 
transformations, a transition from the Higgs phase to a dis- 
tinct confinement phase is possible. The Higgs phase is the 
phase of spontaneously broken center symmetry, while con- 
finement corresponds to the symmetric phase. This division 
of the P — y phase diagram into two separate phases was veri- 
fied numericaUy long ago, in the Monte Carlo investigation of 
ref. [12], which mapped out the approximate location of the 
transition line. 

The Higgs phase of the adjoint Higgs model is often char- 
acterized as a spontaneous breaking of the local gauge sym- 
metry, from SU (2) down to t/(l). In view of the Elitzur theo- 




FIG. 5: Remnant syimnetry order parameter Q vs. P at three values 
of yin the gauge-adjoint Higgs model. 



rem, which states that a local symmetry cannot break sponta- 
neously, this characterization is a little misleading. However, 
as we have discussed above, there exists in Coulomb gauge a 
remnant gauge symmetry which is global on a time shce, and 
which can break spontaneously on a time sUce in the sense 
described in the previous section. We have therefore studied 
the phase diagram of the adjoint Higgs theory via two observ- 
ables: (i) the plaquette energy 



(3.4) 



and (ii) the remnant symmetry breaking order parameter Q 
defined in eq. (2.8). What we find is that the transition lines 
(Fig. 3) detected by each of these two parameters coincide; 
the common line location agrees with the earlier results of 
Brower et al. based on the plaquette energy alone. In Fig. 
4 we plot Ep vs. P at three values of y. The existence of a 
phase transition for the two larger values of y is clearly visible; 
there is no transition apparent at the smallest y value. Fig. 5 
is the corresponding plot of Q vs. P at the same three values 
of y. At P,y values where Ep shows a transition, the transition 
in Q is even more evident. Conversely, where no transition 
is seen in Ep, at y = 0.5, neither is there a transition in Q. 

-1/2 

Finally, in Fig. 6, Q is plotted against V^^ , and we show the 
extrapolation of g to a small value (consistent with zero) at 
infinite volume, for couplings in the confined phase. 



IV. HIGH-TEMPERATURE DECONFINEMENT 

We have seen that in the SU(2)-adjoint Higgs model, things 
go much as one might have expected a priori: remnant sym- 
metry breaking coincides with Z2 center symmetry breaking, 
and in consequence the presence of a confining Coulomb po- 
tential is correlated with the presence of a confining static 
quark potential. One might then guess that renmant and center 
symmetry breaking always go together. This appears not to be 
true, as we have discovered in our investigation of pure SU(2) 
lattice gauge theory at high temperature. 
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SU(2) with adjoint Higgs, confined phase V(R,0)=-log[G(R,1 )], p=2.3 




FIG. 6: Plot of Q vs. root inverse 3-volume, and extrapolation of Q 
to infinite volume, at several couplings in the confined phase of the 
gauge-adjoint Higgs model. 



FIG. S:V{R, 0) in the deconfined phase, at P = 2.3 with n,=2 lattice 
spacings in the time direction, and space volumes ranging from 12^ 
to 323. 



Pure SU(2), |3=2.3, n,=2 
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FIG. 7: The Q parameter vs. root inverse 3-volume in the high- 
temperature deconfined phase, pure SU(2) gauge theory at P = 2.3 
and lit = 2 lattice spacings. 



Monte Carlo simulations of the pure SU(2) gauge theory 
were carried out on x 2 lattices at (3 = 2.3, which is inside 
the deconfined phase. Figure 7 is a plot of Q vs. l/y/Vs, where 
it seems that Q tends to zero at large volume (although we 
cannot entirely rule out a small non-zero intercept at V3 °o). 
This impression is strengthened by our plot of V (/?,0) in Fig. 
8, where it is clear that the Coulomb potential goes asymptot- 
ically to a straight line, as the lattice volume increases. There 
is no indication of the screening of the static quark potential 
(as measured by Polyakov Une correlators), which occurs at 
much smaller distances. 

The results for the Coulomb potential in the deconfined 
phase are not paradoxical; we have already noted that remnant 
symmetry breaking is a necessary but not sufficient condition 
for confinement, and that the Coulomb potential is only an up- 
per bound on the static quark potential. Thus it is possible for 
the Coulomb potential to increase linearly even if the static 
quark potential is screened, as evidently occurs in the decon- 
fined phase. Nevertheless, this result is a little surprising, and 
it would be nice to understand it a little better. 



Recall that in the continuum, the instantaneous part of 
the timelike gluon propagator, which is proportional to the 
Coulomb interaction energy between static charges, is given 
by eq. (2.1 1). Note that this is the expectation value of an 
operator which depends only on the space components of the 
vector potential at a fixed time. On the lattice, this translates 
into an operator which depends only on spacelike links on a 
time-slice. However, we know that spacelike links, at fixed 
time, are a confining ensemble, in the sense that spacelike Wil- 
son loops have an area law falloff even in the high-temperature 
deconfinement phase. If the confining property of the space- 
like Unks on a timeslice is not removed by the deconfinement 
transition, then it is perhaps less surprising that the confining 
property of the (latticized) operator in eq. (2.11), which de- 
pends only on spaceUke Unks on a timeslice, survives in the 
deconfinement regime. 

As a check, we apply a procedure that is known to remove 
the confining properties of lattice configurations. This is the 
de Forcrand-D'Eha [6] method of center vortex removal. The 
procedure is to first fix a given thermalized lattice configu- 
ration to direct maximal center gauge, i.e. the gauge which 
maximizes 

^ = I(2Tr[f/,W])' (4.1) 
and carry out center projection 

Z^(x)=signTr[f/^(x)] (4.2) 

to locate the vortices. Vortices are then "removed" from the 
original configuration by setting 

M^) ^ = Znmix). (4.3) 

In effect this procedure superimposes a thin Z2 vortex inside 
the thick SU(2) center vortices. The effect of the thin vortex 
is to cancel out the long range influence of the thick vortex 
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FIG. 9: The effect of vortex removal on V{R,0) in the deconfined 
phase, P = 2.3 on a 16^ x 2 lattice. 



on Wilson loops. It was found that this procedure not only 
removes the Wilson loop area law falloff, but also removes 
chiral symmetry breaking, and sends every configuration into 
the zero topological charge sector [6]. 

Having removed center vortices from the x 2 lattice, 
thereby removing the area-law falloff of spacelike Wilson 
loops, we fix the modified configuration to Coulomb gauge, 
and compute timelike link correlators in order to measure 
V{R,0). The effect is quite dramatic. It was found in ref. 
[2] that vortex removal in pure gauge theory, in the low- 
temperature confined phase, removes the confinement prop- 
erty of the Coulomb potential. Now we see, from Fig. 9, 
that vortex removal also removes the confining property of 
the Coulomb potential in the high-temperature deconfinement 
phase. This is in accord with the idea that it is the confining 
property of the ensemble of spacelike links at fixed time (or, 
perhaps, the percolation of center vortices on any time slice) 
that is crucial for the confining property of the Coulomb en- 
ergy. 

It is interesting to ask whether there is some connection 
between the center vortex confinement mechanism, and other 
proposals based on the Gribov horizon [14]. While this 
question cannot yet be answered definitively, there is one very 
intriguing fact that may be relevant: when gauge-transformed 
to the minimal Coulomb gauge, thin- vortex configurations lie 
on the Gribov horizon, as we now explain. 



continuum gauge theory, A is convex and bounded in every 
direction [15]. By simple entropy considerations, the popu- 
lation in a bounded region of a high-dimensional space gets 
concentrated at the boundary. For example inside a sphere 
of radius R in a D-dimensional space, the radial density is 
given by r^^ ' dr and, for r<R, is highly concentrated near the 
boundary r = R for D large. With lattice discretization the di- 
mension D of configuration space diverges like the volume V 
of the lattice. Moreover on the boundary the Faddccv-Popov 
operator M(A) = — V-Z)(A) has a vanishing eigenvalue, as we 
will see. (At large volume it has a high density of small pos- 
itive eigenvalues.) This makes the color-Coulomb interaction 
kernel K{x,y;A) = ^ {A){-V^)M^ ' (A) 1;^,. of long range for 
typical configurations A that dominate the functional integral. 
We thus expect thatVcoM/(|x — y|) = {K{x,y;A)) is long range, 
although this qualitative argument is not precise enough to es- 
tablish that Vcoui {R) rises linearly at large R, as suggested by 
the numerical data. 

We shall show that vortex dominance, which is strongly 
supported by the data just presented, is consistent with this 
simple confinement scenario in minimal Coulomb gauge. 
More precisely we shall show that when a center config- 
uration {defined below) is gauge transformed to minimal 
Coulomb gauge it lies on the boundary 3A of the fundamen- 
tal modular region A. According to the confinement sce- 
nario in minimal Coulomb gauge, the probabiUty measure is 
dominated by points at or near the boundary 3A. So center 
dominance, when translated into the minimal Coulomb gauge, 
means dominance by a subset of configurations on the bound- 
ary 3A. This is a stronger condition than the confinement sce- 
nario in minimal Coulomb gauge, but consistent with it. 

Proof of assertion. To simplify the kinematics we give 
the continuum version of the argument. Numerical gauge fix- 
ing to minimal Coulomb gauge corresponds to minimizing on 
each time slice the functional, F{A) = ||A|p, with respect to 
local gauge transformations. Here ||A|p = / d^x |Ap is the 
square Hilbert norm of A'f. At a minimum, which may be rel- 
ative or absolute, the first variation with respect to infinites- 
imal gauge transformations 5A, = £)i(A)a), vanishes for all 
CO, 5||A||2 = 2{Ai,Di{A)03) = 2(A,-,8,co) = 0, which gives the 
Coulomb gauge condition d/A, = 0. Moreover at a relative or 
absolute minimum the second variation with respect to gauge 
transformations 8^ | |A 1 1^ = 2(Z)(A),co, 3;Co) > is non-negative 
for all (0, which is the statement that the Faddeev-Popov oper- 
atorM(A) = — V D(A) is non-negative. These two conditions 
define the Gribov region ii. 



V. CONFINEMENT SCENARIO IN COULOMB GAUGE 
AND VORTEX DOMINANCE 



£1 = {A: diAi = and - 8iA(A) > 0}. (5.1) 

The fundamental modular region A is the set of absolute min- 
ima with respect to gauge transformations. 



In the confined phase [2], and in the deconfined phase. 
Fig. 8, one sees clear evidence of a linearly rising color- 
Coulomb potential. There is a simple intuitive scenario in 
the minimal Coulomb gauge that explains why Vcouiifi) is 
long range [14]. In minimal Coulomb gauge the gauge- 
fixed configurations are (3-dimensionally) transverse config- 
urations that lie in the fundamental modular region A. In 



A={A: ||A|| < ||«A|| foraUg}. 



(5.2) 



It is included in the Gribov region, A C fi. In the interior of 
n all eigenvalues of M(A) are strictly positive X„ > (apart 
from the trivial null eigenvalue associated with constant gauge 
transformations 8,(0 = 0), and on the boundary d£l there is 
a non-trivial null eigenvector 3,Z);(A)coo = 0, and all other 
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eigenvalues are non-negative.^ 

We call a center configuration any lattice configuration 
Z, (x) for which every Unk variable is a center element, Z, (x) e 
Z for every link (x, i). The only non-zero action excitations of 
center configurations are thin center vortices. Such config- 
urations are invariant under all global gauge transformations 
g~^Zi{x)g = Ziix). Now apply an arbitrary local gauge trans- 
formation h{x) to the center configuration Z,(x) V,(x) = 
{x)Zi{x)h{x + 1). We shall take h{x) to be the gauge trans- 
formation that brings the center configuration into the mini- 
mal Coulomb gauge. In general, the transformed configura- 
tion Viix) is not an element of the center, but it is invariant, 
Viix) = g'^ {x)Vi{x)g'{x + r), with respect to the gauge trans- 
formation g' (x) = /z" ' {x)gh{x) which, in general, is no longer 
global. 

We give an infinitesimal characterization of the invariance 
of the configuration Vi{x) under the gauge transformations 
^'(x). The set of global gauge transformations form the SU(A^) 
Lie group and the g'{x) = exp[co(x)] form a representation 
of this group. Here (o(x) is an element of the Lie alge- 
bra of SU(A^). This algebra has — I linearly indepen- 
dent elements co"(x), where n = \,...,N^ — \, that satisfy 
[ro'(x),co'"(x)] = /''™co"(x). Thus the configuration Vi{x) = 
exp[A,(x)], which is the gauge transform of the center con- 
figuration Z, (x) into the minimal Coulomb gauge, is invariant 
under local gauge transformations with — 1 independent 
generators co".^ This is the statement that in continuum nota- 
tion readsA, = Ai + eA (A)(o", or A (A)co" = 0. It follows that 
the CO" also satisfy the weaker condilion V •D(A)(n" = 0. Here 
A is the representative in minimal Coulomb gauge of the cen- 
ter configuration, and as such it lies in the fundamental mod- 
ular region A G A (by definition) that moreover is included in 
the Gribov region A C ii, so we have A e il. However the 
equation V • Z)(A)co" = for A e H means that A lies on its 
boundary dCl. With A C i2 it follows that A lies on the bound- 
ary dA of A, at a point where the two boundaries touch. We 
conclude that the gauge transform of a center configuration 
lies on 3A, as asserted.^ 

The argument just given also applies to abelian configura- 
tions, namely a configuration that lies in an abelian subalgebra 
of the Lie algebra. Such a configuration is invariant under a 
global U(l) gauge transformation. When an abelian configu- 
ration is gauge-transformed to the minimal Coulomb gauge, 
it is mapped into a point where the two boundaries 3A and 



' The Gribov region Q. consists of relative minima that are Gribov copies of 
the set A of absolute minima. Our numerical procedure selects the best 
Gribov copy obtained from eight random gauge copies; there is no known 
method for finding the absolute minimum. On theoretical grounds, one 
expects the sensitivity of oitf results to the choice of Gribov copy be small 
[16]. This expectation is susceptible to numerical investigation, as has been 
done recently for the ghost propagator [17]. 

* The gauge orbit of the center configuration is degenerate and has — \ 
fewer dimensions than a generic gauge orbit. 

^ It should be noted that the numerical procedure that we have used to re- 
move center vortices is 4-dimensional. The argument given here applies 
in Coulomb gauge to center configurations defined within 3-dimensional 
time-slices. 



3i2 touch. But now the equation D,(A)a) = has, in general, 
only one linearly independent (non-trivial) solution, instead 
of A'^ — 1 . Thus abelian dominance is also compatible with 
dominance of configurations on the boundary 8A. 



VI. SU(2) GAUGE-FUNDAMENTAL fflGGS THEORY 

We now consider SU(2) gauge theory with the radially 
frozen Higgs field in the fundamental representation. For the 
SU(2) gauge group, the lattice action can be written in the 
form [18] 

S = P £ iTr[t/t/t/+f/'''] 

plaq 

+ yE5Tr[^^(x)f/^(x)^(x+Ai)] (6.1) 

x,ti 

with (|) an SU(2) group-valued field. This theory cannot be 
truly confining for non-zero y, since the matter field can screen 
any charge, and this simply reflects the absence of a non-trivial 
global center symmetry. However, at sufficiently small y there 
exists a "pseudo-confinement" region, where the static poten- 
tial (as measured by the correlator of Polyakov loops) is linear 
for some intermediate range of quark separations before the 
onset of screening. At large y there is a Higgs region, where 
the linear potential is completely absent. It was shown many 
years ago by Fradkin and Shenker [5] that any two points in 
the Higgs and pseudo-confinement regions can be joined by 
a path in the P — y coupling plane that avoids aU thermody- 
namic singularities. Although there exists a line of first-order 
phase transitions in the P — y plane, this line has an endpoint 
and does not divide the diagram into thermodynamically sep- 
arate phases. (It appears as the soUd Une in Fig. 12, below.) In 
accord with the Fradkin-Shenker observation, numerical sim- 
ulations suggest that there is only one non-confining phase in 
the gauge-fundamental Higgs theory. 

In apparent contradiction to this fact, one can make a strong 
case for the existence of a remnant symmetry-breaking transi- 
tion at small p. We set P = so the action is 

= y£^Tr[^t(^)f;^(^)(^(^+-)]_ (6,2) 

x,ii 

We shall show that (i) at large y there is spontaneous break- 
ing of the remnant gauge symmetry, associated with the short 
range of Vcoui{R), but (ii) at small y there is a linear rise of 
yamiiR) at large R. Thus along the line p = in the P - y 
plane we expect a transition at some finite value of y. 

(i) First consider the limit y — » oo. The action Sa^ is a max- 
imum when f/u(x) = (|)(x)(|)^(x + ju) holds on each link, and 
when y gets large, U^{x) gets frozen at this value. This con- 
figuration is a gauge-transform of the identity U^{x) = /, and 
when fixed to the minimal Coulomb gauge the spatial compo- 
nents get fixed to the identity Ui{x,t) = (|)(x,/)(|)^(x + i,t) = I, 
for all X and i. Thus ^{x,t) is independent of x. We write 
(|)(x,f) = g{t), and we have Uo{x,t) = g{t)g'^{t + 1), which is 
also independent of x. This is the gauge analog of all spins 
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aligned, and we expect spontaneous breaking of the remnant 
gauge symmetry. Indeed, we have 



u{t) 



^i:t/o(x,o^ 



■■g{t)g\t+\), 



(6.3) 



and the order parameter Q defined above has the value 2=1. 
This is maximal breaking of the remnant gauge symmetry. 

(ii) Now consider small values of y. We shall calculate the 
lattice analog of Vcoui{R), namely V{R,0) = -log[G(/?, 1)], 
to leading non-zero order in y. The gauge fixing involves only 
the spatial link variables f/,, and with the action the integra- 
tion over the Uq factorizes into a product over link integrals. 
To evaluate G(R,l) = (iTr[f/o(x, l)t/d(y, 1)]), we first inte- 
grate over the Uo{x,t), with the result, to leading order in y. 



G{R,l) 



t 
16 



Tr 



^t(x,2)^(X,iy(y^2)Ky^l) 



(6.4) 



There are now 4 unsaturated (|) fields. For simplicity we 
suppose that (x, 1) and (y, 1) are joined by a principle axis, 
which we take to be the 1-axis, and R= |x — y|. The lead- 
ing contribution to the (|)-integration at small y is obtained 
by saturating each hnk on the fine that runs from (x, 1) to 
(y, 1) by "bringing down" from the exponent the term 
2yTr[(|)(z,l)f/i(z,l)(|)'''(z + T, 1)], and likewise for the line 
from (x,2) to (y,2). This gives a factor y^^. The ^- 
integrations are now effected. The remaining integration on 
the spatial link variables J/, is finite because of the gauge fix- 
ing. We cannot evaluate it explicitly, but this last integra- 
tion does not introduce any further y-dependence. We thus 
obtain G{R,l) = y^^+2 x H{R). Here H{R) is not known, 
but it is independent of y. This gives V{R,Q) = —{2R + 
2) logy -log[// (/?)]. The asymptotic fall-off in the correla- 
tor G{R, 1 ) should not be more rapid than exponential, so 
lim/f^oo log[H{R)] = —cR. We thus obtain a linear rise at large 
RmV{R,0) GqquiR^ where the "Coulomb" string tension is 
given by 



Ocoh; = -21ogy-Fc. 



(6.5) 



This is non-zero for small y. To make this argument rigor- 
ous one would have to show that the expansion of G(/?, 1) in 
powers of y converges. However this calculation does strongly 
suggest that at P = the remnant symmetry is unbroken for 
small y whereas, as we have seen, it is broken at large y.^ 

Returning to the SU(2) gauge-fundamental Higgs theory, 
we note that the above calculations at large and small y are 
easily understood in terms of the Coulomb-gauge and center- 
dominance scenarios. Indeed, for large y (and any P), the 
gauge- fixed configurations are at or near t/j(x) = /. This con- 
figuration is an interior point of the fundamental modular re- 
gion A. Thus at large y the coupUng to the fundamental Higgs 



' One can make a similar calculation of V(i?,0) for small p (strong coupling) 
in pure SU(A') gauge theory. This yields a finite "Coulomb" string tension 
at small p given by Ocou/ = ~ log P- This suggests that at least in the strong- 
coupling region there is a "Coulomb flux tube" that connects the external 
sources. 
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FIG. 10: Plot of Q vs. root inverse 3-volume, and extrapolation of Q 

to infinite volume in the gauge-fundamental Higgs theory at p = 2. 1 
and various above (y= 1-2) and below the first-order transition 
point around y = 0.9. 



is effective in keeping configurations away from the bound- 
ary 3A, where the thin vortex configurations are to be found, 
and where the Faddeev-Popov operator is of long range. On 
the contrary, at small y the coupling to the Higgs field is inef- 
fective, and entropy leads to dominance of configurations on 
the boundary 3A. 

Now we turn to numerical simulation. Figure 10 is a calcu- 
lation of 2 vs. y at P = 2.1, where it is known (from the work 
of Lang et al. [18]) that the first-order transition is around 
y = 0.9. Below y = 0.9, Q seems to extrapolate to zero at 
infinite volume, while above the transition Q extrapolates to a 
non-zero constant. There appears to be an actual discontinu- 
ity in Q across the whole fine of first-order (thermodynamic) 
phase transitions in the P — y coupling plane. At sufficiently 
small values of P, there is no thermodynamic transition, and 
we see no discontinuity in 2 as a function of y. What we see 
instead is that 2 « over a finite range of y, and then, beyond 
a critical value y = y^r, Q smoothly increases with increasing 
y. Our results for g vs. y at p = 0, on 8"* and 16"* lattices, are 
shown in Fig. 1 1 ; the solid line is the presumed extrapolation 
to infinite volume. If we were dealing with a spin system, and 
Q were the magnetization, this would clearly represent a sec- 
ond order phase transition. In the present case it is certainly a 
synmietry-breaking transition, separating a symmetric region 
with 2 = from a broken- symmetry region of 2 > 0- 

On the other hand, despite the existence of a symmetry- 
breaking transition, there is no thermodynamic transition of 
any kind at P = 0. At this value of P the free energy can be 
computed exactly, with the result, in a 4- volume V 



F(y) =4yiog 



2/1 (y) 



(6.6) 



which is perfectly analytic for all y > 0. Thus we have con- 
firmed the theoretical argument that there must be a remnant 
symmetry-breaking transition even at small P, but we have 
also found that this transition is not accompanied (at small 
P) by a thermodynamic transition, defined as some degree of 
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FIG. II: 2 vs. y at P = in the SU(2) fundamental Higgs model, on 
8^ and 16^ lattices. The solid line is the presumed extrapolation of Q 
to infinite volume. 
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FIG. 12: Phase diagram of the SU(2) fundamental Higgs 
model. There is a thermodynamic transition and a Q (remnant 
symmetry-breaking) transition along the solid line, but only a non- 
thermodynamic transition (Kertesz line) along the dashed line. 



non-analyticity in the free energy. Our result for the Une of 
critical couplings of the remnant symmetry-breaking transi- 
tion is shown in Fig. 12. Along the solid line there is also a 
thermodynamic (first-order) transition, which is absent along 
the dashed line. 

Our result for the gauge fundamental Higgs system is not 
entirely new; it was in fact anticipated by Langfeld in ref. [3], 
who considered a closely related model in Landau (rather than 
Coulomb) gauge. In that work the modulus of the space-time 
averaged Higgs field was used as an order parameter to detect 
the breaking of remnant symmetry, which in Landau gauge 
must be both space and time independent. In ref. [4] it was 
further suggested that the line of remnant symmetry breaking 
transitions, where it is unaccompanied by a line of thermo- 
dynamic transitions, is a Kertesz line [19]. A Kertesz line is 
a Une of percolation transitions; the original example comes 
from the Ising model. In the Ising model, in the absence of 
an external magnetic field, there is a phase transition from a 
Z2 symmetric to an ordered phase, and this transition can be 



expressed, in different variables, as a transition from a perco- 
lating phase at low temperature, to a non-percolating phase at 
high temperature. In the presence of a magnetic field, the par- 
tition function and thermodynamic observables become ana- 
lytic in temperature; there is no thermodynamic phase tran- 
sition. Nevertheless, the percolation transition persists, and 
traces out a Kertesz Une in the temperature-magnetic field 
plane, completely separating the phase diagram into two re- 
gions. But if there is a Kertesz line in the gauge-Higgs cou- 
pling plane, the question is what kind of objects are percolat- 
ing. Based in part on results reported by Bertie and Faber [20], 
Langfeld [4] proposed that the unbroken remnant symmetry 
region is a region of percolating center vortices, which cease 
percolating in the broken symmetry region. There is now very 
good evidence for a vortex percolation transition of this kind 
in gauge-fundamental Higgs theory, reported in ref. [21].^ 

Our findings here support the idea that there is some phys- 
ical distinction that can be made between the Higgs and 
the pseudo-confining regions of the gauge-fundamental Higgs 
phase diagram. In the pseudo-confining region the remnant 
symmetry is unbroken, the Coulomb potential rises Unearly, 
and center vortices percolate, while the Higgs region is a re- 
gion of broken symmetry, the Coulomb potential is asymptoti- 
cally flat, and center vortices do not percolate. This distinction 
appears to exist despite the fact that the two regions are ther- 
modynamically connected, as demonstrated by Fradkin and 
Shenker in ref. [5]. 

Before leaving the gauge-fundamental Higgs theory, we 
should ask what are the effects of vortex removal in the sym- 
metric phase. We already know that the Coulomb potential 
in the fundamental Higgs theory must be confining in the 
syimnetric (pseudo-confinement) phase and screened in the 
broken phase; the only issue is how the Coulomb potential 
is affected in each phase if center vortices are removed. In 
Figs. 13 and 14 we see the Coulomb potential in the symmet- 
ric (P = 2.1, y = 0.6) and broken (P = 2.1, y = 1.2) phases, 
respectively, before and after vortex removal. In the symmet- 
ric phase, vortex removal by the de Forcrand-D'Elia proce- 
dure sends the Coulomb string tension to zero, as in the high- 
temperature phase of pure gauge theory. Deep in the Higgs 
phase, on the other hand, the effect of vortex removal is seen 
to be very minor. 



VII. THE ADJOINT-HIGGS MODEL REVISITED 

We have understood theoretically why there should be a 
remnant symmetry-breaking transition at small P in the gauge- 
fundamental Higgs model, as seen in the numerical data. This 
raises immediately the question of why there is no correspond- 
ing transition in Q, at small P, in the gauge-adjoint Higgs the- 
ory. Now in the fundamental-Higgs model, the transition at 
small P appears to be a percolation transition, as discussed in 



' The role of Kertesz Unes in high-temperature QCD is also discussed by 
Satz in ref. [22]. 
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FIG. 13: Gauge-fundamental Higgs theory: effect of vortex removal 
in the symmetric (pseudo-confinement) phase, P = 2. 1 , y = 0.6. 
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phase, in gauge-adjoint Higgs theory. 



the previous section. But the absence of a percolation transi- 
tion at low P in the adjoint-Higgs case is easy to understand. 
There the Higgs part of the action is insensitive to the exis- 
tence of thin center vortices, and cannot suppress their con- 
densation at small P regardless of the value of j. More con- 
cretely, at P = 0, the action is invariant with respect to local 
transformations Uo{x) z{x)Uo{x) where z(x) = ±1, and this 
immediately implies that 2 = at infinite volume, again re- 
gardless of the value of y. 

On the other hand, if the absence of a remnant symmetry- 
breaking transition at low P is due to large fluctuations of cen- 
ter elements z(x), then one might still expect breaking, at low 
P and large y, of the SO(3)=SU(2)/Z2 part of the remnant sym- 
metry group, which is insensitive to Uq{x) z{x)Uo{x) fluc- 
tuations. The relevant order parameter is 



Qadj 




(7.1) 



where Uadj{t) is the spatial average of timelike links in the 



adjoint representation 



Uad}{t) = —'£Uo,adj{li,t), 



J-Tr 



a"Uo{x,t)a''U^{x,t) 



(7.2) 



Rather surprisingly, there appears to be no transition in Qadj 
either, in the adjoint-Higgs model. As we see in Fig. 15, for 
data taken on a small S'* lattice, there is no sign of any tran- 
sition for Qadj at flnite y and P = 0. In fact, Qadj is essen- 
tially y-independent. Extrapolation of Qadj to infinite volume 
anywhere in the confined (Q = 0) phase is consistent with 
Qadj = 0, as seen in Fig. 16. Thus the transition line for Qadj 
in the P — y coupling plane appears to be the same as for Q 
in Fig. 3, but not like Q in Fig. 12, and this fact calls for an 
explanation. 

The difference in the phase diagrams of Figs. 3 and 12 oc- 
curs at small P and large y. To understand this difference we 



13 



set P = in the adjoint Higgs action (3.2) so it reads 



5* = |lrW<./(^+Ai), 



(7.3) 



and we shall evaluate the equal time correlator in minimal 
Coulomb gauge 



(7.4) 



for Y large. 

Step 1. Integration over Uo- The gauge fixing does not 
involve the Uq, so the integrals over the Uo factorize, with the 
result on each time-Uke Unk 



^ ^ cosh(|) 2 
sinh(l) y' 



(7.5) 



Here X has the limiting value A, = 1 at large Y- The correla- 
tor factorizes into a product of expectation values on adjacent 
time slices 



Gadj{x-y, 1) = i (r (x, l)f (x,2)f' (y,2)r (y, 1) 



= i(r(x,l)r(y,l)) (^^x,2)^^(y,2) 
This gives 

Gadjix-y,\) = 3C\x-y), 

where 



C(x-y) = i(r(x,l)r(y,l)) 



(7.6) 



(7.7) 



(7.8) 



is the (|)-propagator within a single time sUce. It is calculated 
using the action within a time slice, 



Sslice = |£r(x)f/^*^/x)^^(x + f). (7.9) 



(The argument / that enumerates time slices is suppressed in 
the following.) 

Step 2. Reduction to U(l) theory. To evaluate the (|)- 
propagator at large Y^ we introduce a gauge transforma- 
tion ^(x) e SU(2) that depends on (|)(x) with the property 
that the group element in the adjoint representation, denoted 
Hl^'^lsix)], is the SO(3) rotation matiix that rotates the 3- 
direction into the (|)-direction, ^"{x) = ^^''•'[^(x)]. The action 
and correlator are given by 



Sslice = ]:Y,'J{.'' [g-\x)Ui{x)g{x + l) 

C(x-y) = i(^«3[^(x)]^«3[g(y)]). 



(7.10) 



These expressions are invariant under g{x) ^ g{x)h{x), where 
h{x) is an SU(2) element that corresponds to a rotation 
about the 3-axis. So we may average over h{x) which re- 
sults in replacing the integral over (|)(x) by an integral over 



g{x) E SU(2). At large Y the functional integral is domi- 
nated by the maximum of the action SsUce, which occurs where 
li^^[g~^{x)Ui{x)g{x + i)] = 1 holds on each link, namely, 
where Ui{x) = g~^(x)Ui{x)g{x -|- i) is a rotation about the 3- 
axis. These rotations form the U(l) group. Thus at the maxi- 
mum of the action the link variables are given by 



Ui{x) = g{x)uiix)g \x + i) 

= ^Ui{x), 



(7.11) 



where g(x) e SU(2). Thus f/,(x) is an SU(2) gauge transform 
of M,(x) e U(l). At large Ythe Ui{x) get frozen into this form, 
and the integral over Ui{x) G SU(2) gets reduced to an integral 
over abeUan configurations m,(x) e U(l). We have noted that 
we may replace g{x) by g{x)h{x), where h{x) E U(l). We use 
this freedom to gauge fix the Mj(x) within the U(l) group of 
rotations about the 3-axis. Naturally we choose the minimal 
Coulomb gauge of U(l) gauge theory. Thus the integral over 
the Ui{x) e SU(2) gets replaced by an integral over gauge- 
fixed configurations m;(x) S U(1). We designate this set by 
T. In continuum gauge theory, T is the set of all transverse 
abelian configurations A^(x). In sharp contrast to the SU(2) 
case, in abelian gauge theory there are no Gribov copies. Dif- 
ferent transverse configurations A? (x) are gauge-inequivalent, 
and the set T of gauge-fixed abeUan configurations A? (x) is 
the unbounded set of all transverse abelian configurations, 
without restriction to a fundamental modular region. (This 
is a major difference between abelian and non-abelian gauge 
theories, and is the basis of the confinement scenario in mini- 
mal Coulomb gauge.) 

Step 3. Integration over the g{x). Although u G T is 
completely gauge-fixed by minimizing with respect to local 
t/(l) gauge transformations, it is not gauge-fixed by min- 
imizing with respect to local SU(2) gauge transformations, 
because that is a larger group. But our calculation in min- 
imal Coulomb gauge requires that the SU(2) configuration 
Ui{x) = *M,(x) e A be completely gauge-fixed with respect 
to local SU(2) gauge transformations. So for given u, g{x) is 
the unique SU(2) gauge transformation (modulo global gauge 
transformations) that accomplishes this gauge fixing. We 
write g{x) = g{x;u), and the (|)-propagator is given by 

C{x-y) = j^duP{u) \%,"Mx-u)]'K,''\g{y;u)], (7.12) 

where P{u) is a positive probability density. 

Step 4. Integration over the u{x). We saw in section V 
that when configurations in an abelian subgroup are gauge- 
fixed, they get mapped into the boundary 8A. However that 
conclusion is an over- statement which ignores the fact that 
some U(l) configurations u E T lie in the interior of A (al- 
though this happens with probabihty 0, as we shall see). Since 
in the continuum Umit A is a subset of transverse SU(2) con- 
figurations A"'"'{x) that is bounded in every direction whereas 
T is the unbounded set of all transverse abelian configurations 
A^'"^{x), it follows that some configurations uGT lie inside A 
and some lie outside A. Correspondingly we break up the in- 
tegral into contributions from uET inside A, and from u eT 
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outside A, 

C(x-y) 
Qn(x-y) 

CoM!(x-y) 



Q„(x-y)+Co„;(x-y) 

/ jMP(M)i3l''3[^(x;M)]3l«3[^(y.„)] 

JrnA 



du P{u) X 

T-TnA 



(7.13) 



For u E T OA the unique gauge transformation g{x,u) that 
brings u inside A is the identity transformation 

g{x,u)=r, uGTnA. (7.14) 

This gives ^''^[g{x;u)] = = 8a3, and we obtain 

Q„(x-y) = / duP{u) 
JuernA 

= Pin, (7.15) 

where Pi„ is that probability that u GT lies inside A. This is 
independent of x and y, and corresponds to ordered spins. It 
resembles the calculation with coupling to fundamental Higgs 
at P = and y large, where Pin = 1. For u outside A it ap- 
pears that the solution ^(x,m) of the spin-glass minimization 
problem depends in a very irregular and disordered way on u 
and X, so Cout{R) decays rapidly, Um^^oo Coaj (i?) = 0. With 
GadjiR, 1) = 3C2(/?) this gives 



^mGadj{R,\) = 2>Pl 



(7.16) 



Thus, remarkably, a numerical determination of GadjiRA) 
provides a direct measurement of the probability Pi„ that a 
configuration u eT lies inside A. 

The data of Fig. 16 strongly suggest that Qadj extrapolates 
to at infinite volume. This is the disordered phase, in which 
lim/;_„,o Gadj{R, 1) = holds. This gives 



Jim P;, 



:0. 



(7.17) 



We have noted that in continuum gauge theory A is bounded 
in every direction whereas T is unbounded in all directions. 
According to the simple entropy estimate, p{r)dr ^ r^^^dr, 
where D is the (very high) dimension of the space T of trans- 
verse U(l) configurations, the fraction P;„ of the set T that 
lies inside the bounded region A is negligible, f,„ — > 0, as 
in fact the data indicate, and we have also Pout 1 • More- 
over an abelian configuration u E T that lies outside A, gets 
gauge-transformed in the minimal Coulomb gauge into a con- 
figuration U = ^u that lies on the boundary 3 A, as was shown 
in section V, so in this instance all the probability lies on the 
boundary 3A and the measure of the interior vanishes. This 
exemplifies the simple scenario in Coulomb gauge, according 
to which confinement occurs when the functional integral is 
dominated by the boundary 3A. 

The absence at P = of a transition in Qadj as y increases 
from to oo is now explained. For with coupling to the ad- 
joint Higgs, the measure of ordered spins that would break 
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FIG. 17: Plot of Q vs. root inverse 3-volume, and extrapolation of 
Q to infinite volume in QED4, for P = 0.7 (confining phase) and 
P = 1.3 (massless phase). 



the remnant gauge symmetry is Pin = at large y, and the rem- 
nant gauge symmetry is preserved. By contrast, with coupling 
to the fundamental Higgs for y large, Ui{x) is a gauge trans- 
form of the identity /, as we have seen in section VI. Since 
/ is certainly in A, then P,„ = 1, and the remnant symmetry 
is maximally broken. The coupUng to the fundamental Higgs 
field at large y keeps configurations away from the boundary 

aA. 

Finally we wish to emphasize that with couphng to the ad- 
joint Higgs, the numerical result Pi„ = is a direct manifes- 
tation of the deep difference between the fundamental modu- 
lar regions (defined by minimizing F[A] with respect to local 
gauge transformations) of an abelian and non-abelian gauge 
theory. The fundamental modular region is unbounded in ev- 
ery direction in an abelian gauge theory, but bounded in every 
direction in a non-abelian gauge theory. 



Vni. COMPACT QED 

Although in the preceding sections we have considered only 
the theories with a non-abelian SU(2) gauge invariance, there 
is no barrier to studying remnant symmetry breaking in an 
abelian model such as compact QED 4. In this model we have 
confinement at strong couplings and a massless phase at weak 
couplings, with a transition between the two phases at approx- 
imately P = 1. In Fig. 17 we show our results for (the abelian 
analogue of) 2 at P = 0.7, which is inside the confinement 
phase, and at P = 1.3, which lies in the massless phase. The 
results are as expected; Q extrapolates nicely to zero at infi- 
nite volume in the confined phase, and appears to extrapolate 
to a non-zero value in the massless phase. Figure 18 is a plot 
of V(/?,0) vs. charge separation /? at P = 0.7 on variety of hy- 
percubic lattice volumes L^. Note that the potential is insensi- 
tive to changes in lattice volume, and that deviations from the 
linear potential, where they are statistically significant, arise 
from the lack of rotation invariance at strong couplings. 

The situation changes drastically in the massless phase at 
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FIG. 18: The potential V{R,0) in at P = 0.7, on a variety of 

lattice volumes. 



IX. TRANSLATION TO TEMPORAL GAUGE AND 
STRINGLESS STATES 

In this section we shall show how to translate back and forth 
between the minimal Coulomb gauge and the temporal gauge, 
Aq = 0, so the measurements reported here have an equiva- 
lent description in the temporal gauge. Both of these gauges 
are compatible with a Hamiltonian formulation and a physi- 
cal transfer matrix. Moreover the temporal gauge is invariant 
under all space-dependent but time-independent gauge trans- 
formations g{x), where x is a 3-vector, and it may be helpful 
to express things in a more gauge-invariant way. We shall 
see that the state obtained here by numerically gauge-fixing 
to the minimal Coulomb gauge becomes, after translation into 
temporal gauge, a "stringless" state of the type introduced by 
Lavelle and McMuUan [7] that does not involve a Wilson line 
joining the sources. 
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FIG. 19: The potential V{R,0) in QED4 at p = 1.3, on several L'^ 
lattice volumes. 



A. Temporal gauge 

In temporal gauge the continuum Hamiltonian has the 
canonical form 



H=^,jd\{E^+B\ 



(9.1) 



Here the spatial components A"(x) of the vector potential 
are canonically conjugate to the color-electric field operator 

W = 'SSfW ' and Bt = ZijkidjAl + ir^^AjA^). 

The color-charge density of the static quarks does not ap- 
pear in H because Aq = 0. Instead it appears in the generator 
of local space-dependent gauge transformations. 



G«(x) = -Z3-(A).£'^(x) + p^(x), 



(9.2) 



P = 1 .3, where our results for V (/?, 0) are displayed for 5^, 10^ 
and 20"^ lattices in Fig. 19. Although the string tension ex- 
tracted from y(/?,0) appears to be non-zero on each of these 
lattices (which cannot be correct for the Coulomb potential 
in the weak-coupling regime), this string tension drops very 
markedly (by about 1/3) with each doubUng of the lattice vol- 
ume. This volume dependence of the string tension is very 
different from what we have observed for the non-abehan the- 
ories in various phases. For the abelian theory, the expected 
result Of „„/ = is presumably recovered in the large volume 
limit. In fact, from the value of Q extrapolated to infinite vol- 
ume, we can estimate that on a large lattice V (/?,0) should be 
bounded, very roughly, by — ln((2^) « 1.6 at large R. This 
bound impUes that V(^,0) has to level out in large volumes, 
resulting in (Seoul = 0. 

What we learn from the weak-coupUng case is that V{R,Q) 
may have very significant lattice size dependence in theories 
such as QED4 in the massless phase, where the correlation 
length is large (comparable to lattice size) or infinite. It is 
therefore important to compute V{R, 0) on a variety of lattice 
sizes, and to extrapolate Q to infinite volumes, as we have 
done in the preceding sections. 



where 



P"{x) =t"5{x-Xl)+t25{x-X2) 



(9.3) 



is the color-charge density of a pair of external quarks at xi 
and X2- The wave functional ^ap(A) bears the color indices 
of the external quarks, on which the quark charge vectors act 
according to 

(fl ^)ap = fay^Tp; (f2*)ap = -t^;'f'ay, (9.4) 

where the t" form the basis of an irreducible representation 
of the Lie algebra of SV(N), [f",f*] = if'^t". One may ver- 
ify that G(a)) — J d^x (Sf{x)G"-{x) generates an infinitesimal 
gauge transformation. 



[G(co),Af] 
[G(co),£f] 
[G(co),p^] 



= if 



■ahc Tphf^c. 



(9.5) 



This transformation leaves H invariant, [G(co),//] = 0, so we 
may choose wave-functionals that transform irreducibly under 
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the local gauge group. Physical wave-functionals are required 
to satisfy the subsidiary condition 

G(x)>P = 0, (9.6) 

which is both Gauss' law and the statement that the wave- 
functional is gauge-invariant. This condition determines the 
gauge-transformation properties of the wave-functional 

^{'A)=g{xi)^{A)g\x2), (9.7) 

where we use matrix notation for the quark color indices, and 

B. From temporal gauge to minimal Coulomb gauge, and back 

The continuum temporal gauge does not really provide a 
well-defined quantum mechanics because the irmer product 
for gauge-invariant wave functionals 

i^u'Vi) = jdA ^<\{AYV2{A) (9.8) 

diverges due to the local gauge invariance of the wave- 
functionals. Gauge-fixing is required to correctly normalize 
the wave-functionals in temporal gauge. We do this by ap- 
plying the Faddeev-Popov formula to gauge-invariant inner 
products. For this purpose we parametrize configurations by 
A = '^A*, where A^ is the representative of A in the minimal 
Coulomb gauge, so A'' e A is a transverse configuration in the 
fundamental modular region. The Faddeev-Popov formula 
gives 

(^i,*?!) = / dA^ detM(A'') %(A^W2{A^). (9.9) 
Ja 

Here M(A^) = —V • Z)(A"^) is the Faddeev-Popov operator, 
which is symmetric and positive for A"^ e A. The right hand 
side is the iimer product in minimal Coulomb gauge. Thus, in 
the 3-dimensional operator formalism, the minimal Coulomb 
gauge is a gauge-fixing within the temporal gauge of the 3- 
dimensional local gauge invariance, and the wave functional 
in minimal Coulomb gauge is the restriction of the wave- 
functional in temporal gauge to the fundamental modular re- 
gion^^ 

'Pco„;(A'^)='P(A'^); A'^eA. (9.10) 

Conversely, the gauge invariance of the wave-functional 
in temporal gauge defines the unique extension of the wave- 
functional in minimal Coulomb gauge into a wave-functional 
in temporal gauge. We parametrize an arbitrary configuration 



To avoid confusion we note that the two conditions Aq = and 3,A,- = 
can be imposed at a fixed time, which is sufficient for the 3-dimensional 
operator formalism. They cannot both be maintained for all time, and in 
the 4-dimensional Feynman path-integral formalism in Coulomb gauge one 
maintains djA,- = 0, but Aq ^ 0. 



by A = ^A"', where A"' G A, and g{x\A) and A^{x\A) depend 
on the configuration A. The existence and uniqueness of these 
quantities at the non-perturbative level is assured (with lattice 
regularization) by the existence of an absolute minimum with 
respect to gauge transformations of ^a(^) = ||^' A|p. From 
(9.7) above we obtain, in matrix notation, 

'P(A) = 'PPA"'] 

= g{xi;A)'V[A''{A)]g\x2;A), (9.11) 

which expresses the wave-functional in temporal gauge in 
terms of the wave-functional 'P(A'^) in minimal Coulomb 
gauge. The gauge transformation g{x;A) that is found numer- 
ically when gauge-fixing to the minimal Coulomb gauge has 
reappeared in the wave -functional in the temporal gauge. 

[For completeness, we note that a quick way to obtain the 
Hamiltonian in Coulomb gauge [8] is to apply the Faddeev- 
Popov formula to the matrix elements of E^, 

(£f'Pi,£f'P2)= / JA'MetM(A*)(£;'Pi)*£,'P2U=A.r- 

(9.12) 

From this formula the matrix elements of the Hamiltonian in 
Coulomb gauge are easily found once £i'P|^=^tr is specified. 
To evaluate £','P|^^^tr we solve Gauss' law for the longitudinal 
part of the color-electric field. We write Ei Ef — where 
(^^{x) is the color-Coulomb potential operator, and Ef = i^ 
satisfies 

[Ef{x)Af{y)] = i[bij - did0)-'] d{x-y) 5«^ (9.13) 

Gauss' law in temporal gauge reads 

-Didi<^W +Ai X Ef^ = pq^, (9.14) 

where {X x Yf = f^bcxbyc is the Lie bracket. We solve for 
the color-Coulomb potential 

(^{x) ^U=^tr = [M-\A^)p]{x) ^{A^), (9.15) 

where A^ e A. Here M{A^) = -D{A^)-d = -d- D{A^) is the 
Faddeev-Popov operator, and 

p = -Afx£f + p^ (9.16) 

is the total color-charge density of quarks and dynamical glu- 
ons. This gives the desired expression, 

Ei 'PU^^tr = [Ef - diM-\A^)p] ^'(A^), (9.17) 

which is to be used in (9.12).] 



C. Energy in Coulomb gauge is energy in temporal gauge 

The quantity we have measured is the expectation value 

{H,oul) -Eo = Vcoul{\x - y|) + Ese (9.18) 
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in the state with wave-functional in minimal Coulomb gauge 
'I'ap(A^) = 75 5„p3.o(A''), (9.19) 

where 3>o(A'') is the vacuum state of pure glue. 

We translate this back into temporal gauge. From (9. 1 1) we 
obtain 

^ap(^) = fe(xi;A)gt(x2;A)]„p 4>o(A), (9.20) 

where we have used the gauge invariance of the vacuum state 
of pure glue, «I>o(A"') = 4>o(*A"^) = 4>o(A). Thus the quantity 
we measure may equivalently be described as the expectation- 
value of the Hamiltonian in temporal gauge, with this wave 
functional 

('P,//'P) -£o = Vcoui{\^-y\)+Ese. (9.21) 



D. Stringless states 

A wave-functional of type (9.20), with 4>o(A) an arbitrary 
gauge-invariant scalar function, has been considered before 
[7] and was called a "stringless" state. The motivation 
for constructing this state was that it has the correct gauge- 
transformation properties, which is equivalent to Gauss' law 
being satisfied exactly. Moreover in QED, with external 
charges only, the stringless state is the exact wave-functional, 
with g{x\A) = exp[;(V^)^' V • A]. This leads us to expect that 
in QCD the stringless state becomes exact at short distance 
|x — y I . Thus the stringless state has several attractive proper- 
ties. 

The stringless state was originally constructed using a per- 
turbative expansion for g{x;A). This expansion does not con- 
verge for A*^ outside the fundamental modular region A. How- 
ever, as we have noted, the existence of g{x\A) at the non- 
perturbative level is assured by the minimization procedure, 
and g{x;A) is the gauge transformation that we have found 
numerically. 

Since our numerical data strongly suggest that Vcoui(R) 
rises linearly, it appears from (9.20) that the "stringless" state 
of quarks manifests a finite string tension Ocoui at large separa- 
tion. Our numerical finding is that this string tension exceeds 
the standard string tension o by <5coui ~ 3a. This provides a 
measure of the extent to which the stringless state (9.20) fails 
to be the exact ground state of a pair of external quarks. 



It was called "stringless" because in its constraction the thin string, 
Pexp{J^^Ajdx'), was replaced by g{xi;A)g^ {x2;A) which transforms in 
the same way under gauge transformation. It was argued that the thin 
string has infinite energy, whereas the "stringless" state has finite energy 
(after ultraviolet renormalization) and is a better approximation to the cor- 
rect hadron state. However when A lies on a degenerate orbit, ''A = A for 
some h{x) ^ /, the parametrization A = ^A" is singular. As a result the 
stringless wave-functional is singular for such configurations, which may 
raise its energy significantly. 



X. CONCLUSIONS 

In this article we have shown that the confining property of 
the color-Coulomb potential is tied to the unbroken realization 
of a remnant global gauge symmetry in Coulomb gauge. We 
have studied this type of confinement numerically in SU(2) 
gauge-Higgs theories, and in pure gauge theory at zero and at 
finite temperatures. Confinement in the color-Coulomb poten- 
tial is not identical to confinement in the static quark potential. 
We have seen that the deconfined phase in pure gauge the- 
ory, and the pseudo-confinement region of gauge-fundamental 
Higgs theory, are instances in which the color-Coulomb po- 
tential is asymptotically linear, even though the static quark 
potential is screened. In terms of symmetries, the point is that 
center symmetry breaking, spontaneous or expUcit, does not 
necessarily imply remnant symmetry breaking. 

The existence of a confining color-Coulomb potential, in 
cases where the static quark potential is screened, has some 
bearing on the question: In what sense does confinement exist 
in real QCD, with dynamical quarks? The problem is that in 
gauge theories with matter fields in the fundamental represen- 
tation, such as real QCD, there is no non-trivial center symme- 
try, and no possibiUty of having an asymptotically confining 
static potential. Further, in gauge theories with a scalar matter 
field in the fundamental representation, there is no local order 
parameter that can distinguish between the Higgs and confine- 
ment phases, and the Fradkin-Shenker theorem assures us that 
there is no thermodynamic transition of any kind that can iso- 
late the Higgs phase from a distinct confinement phase. AH 
this suggests that there is no fundamental difference between 
real QCD and SU(3) gauge theory in a Higgs phase, and it 
should be possible to interpolate smoothly from the one the- 
ory to the other by a continuous change of parameters. 

The results reported in this paper suggest otherwise. In 
gauge theories with matter fields in the fundamental represen- 
tation, the "confinement" phase and the Higgs phase are dis- 
tinguished by the symmetric or broken realization of a rem- 
nant gauge symmetry. Remnant symmetry breaking is not 
accompanied by non-analyticity in the free energy, neverthe- 
less the (non-local) order parameter Q is directly related to 
the color-Coulomb potential, which is confining in the 2 = 
phase. Thus, in real QCD, the gluon propagator in Coulomb 
gauge is confining. This confining property of (dressed) one- 
gluon exchange is absent in the Higgs phase of a gauge field 
theory. 

We have also uncovered a strong correlation between the 
presence of center vortices and the existence of a confining 
Coulomb potential. Thin center vortices, as pointed out above, 
lie on the Gribov horizon. In every case where we find a con- 
fining Coulomb potential, we also find that removal of center 
vortices takes the Coulomb string tension to zero. In related 



While it is sometimes suggested that confinement should simply be under- 
stood as the condition that asymptotic particle states are all color singlets, 
this condition is also fulfilled in the Higgs phase of gauge-Higgs theories 
(cf., e.g., ref [13]). 
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work on the SU(2) gauge-fundamental Higgs theory, it was 
suggested [4], and recently verified [21], that vortices perco- 
late throughout the lattice in the pseudo-confined phase, and 
do not percolate in the Higgs phase. These phases correspond 
to regions of unbroken and broken remnant symmetry, which 
are completely separated by a percolation transition located 
along a Kertesz line. 

There are many open questions. Presumably the linear 
color-Coulomb potential is associated with a flux tube of lon- 
gitudinal color electric field. If this is really so, does the 
tube have string-like properties; i.e. roughening and a Liischer 
term? Since o„)ui is several times greater than o, the Coulomb 
flux tube must be an excited state. By what mechanism is 
the string tension o of the minimal energy flux tube reduced 
below the value Gcoui of the Coulomb flux tube? Does this 
mechanism involve production of constituent gluons, along 



the lines of the gluon-chain model [23], or is some other pro- 
cess at work? We hope to address some of these issues in a 
future investigation. 
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